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SUMMARY 

The Large Scale Advanced Prop-Fan hub assembly forms a semi-rigid link be- 
tween the blades, which provide the thrust, and the engine shaft, which pro- 
vides the torque. The hub and tailshaft is a one piece partially forged part 
which is carburized, heat treated and machined. A single row ball bearing 
restrains each of the eight blades in the hub, while the tailshaft secures 
the propeller to the engine shaft with two cone seats that are pre loaded 
against each other by the Prop-Fan retaining nut. The hub also forms the 
support for the pitch change actuator system, the control and the spinner. 

The retention transmits the loads from the blades to the hub while allowing 
the changes in blade pitch. The single row ball bearing retention provides 
ease of maintenance by allowing individual blade replacement without disas- 
sembly of the hub. It has a through hardened inner race which seats against 
the aluminum blade shank and an outer race which is integral with the bar- 
rel. The outer race area is carburized to achieve the hardness necessary to 
support the ball loads. The balls are kept from contact with each other by a 
separator. The rotational speed of the propeller keeps the retention sub- 
merged in the oil which is contained in the hub by a seal. 

The hub assembly, because it is the foundation of the blade, also is an inte- 
gral part of the mechanism which determines the natural frequencies of the 
propeller system. Analysis assured that none of these modes intersect with 
the operating frequencies of the Prop-Fan. 

The hub portion of the assembly was designed to the same requirements which 
would be established for a production Prop-Fan (lightweight). Because the 
LAP is destined for test on an existing gearbox, and must accommodate an ex- 
isting propeller control, the tailshaft portion of the assembly is identical 
to the 54H60 propeller tailshaft. Production Prop-Fans could incorporate 
either a tailshaft or a flange dependent on the reduction gearbox configur- 
ation . 
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In recent years, considerable attention has been directed toward improving 
aircraft fuel consumption. Studies have shown that the inherent efficiency 
advantage that turboprop propulsion systems have demonstrated at lower cruise 
speeds may not be extended to the higher speed of today's turbofan and 
turbojet-powered aircraft. To achieve this goal, new propeller designs which 
feature more blades with thin airfoils and aerodynamic sweep are required. 

Since 1975, Hamilton Standard has been deeply involved with the NASA Lewis 
Research Center in the development of the advanced turboprop or Prop-Fan. 

Many aircraft system studies have been accomplished for a variety of subsonic 
air transport applications and all these studies have shown significant fuel 
savings with Prop-Fan propulsion. The fuel savings potential of future 
Prop-Fan powered aircraft is generally 15-20% for commercial applications and 
25-35% for military patrol aircraft compared to equal technology turbofan 
systems, depending upon the specific application, cruise speed, stage length 
and other requirements. 

To date, several models have been designed, manufactured and subjected to a 
number of tests. A series of small-scale 0.6223 meter (24.5 inch) diameter 
module tests have been conducted in both UTRC and NASA wind tunnels and on a 
modified NASA airplane. These tests have shown that propellers with 8-10 
swept blades, high tip speeds and high power loadings can offer increased 
fuel efficiencies at speeds up to 0.8 Mn. 

HSD has designed a 2.743 meter (9-foot) diameter single-rotation Prop-Fan. 
Following the manufacture of the Prop-Fan system, it will be tested at Wright 
Field and in the ONERA S-l wind tunnel in France. The hardware will then be 
used in a follow-on program where it will be run with an engine on a static 
test stand and on a research aircraft. The major objective of this testing 
is to establish the structural integrity of large-scale Prop-Fans of advanced 
construction in addition to the evaluation of aerodynamic performance and 
aero-acoustic design. 

The report which follows covers the design analysis of the Large-Scale Ad- 
vanced Prop-Fan Hub/blade retention assembly. Specifically, analysis of the 
retention area of the blade shank, blade retention, hub, and tail shaft are 
covered. Subjects included are stress and strain analysis, material hardness 
requirements, weight predictions and stiffness characteristics. 
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RETENTION AND HUB LOADS 


UNITED 

TECHNOLOGIES 


A cross section through the hub and retention is shown in Figure 1. All the 
external loads on the retention and hub have their source in the blade 
loads. An understanding of the hub loads requires an explanation of those 
blade loads that are transferred to the hub. The loads on the blade are de- 
rived from the aerodynamic and inertia characteristics of the blade. 

The aerodynamic loads generated by the moving airfoil blade sections are the 
lift and drag shear forces. On a particular blade (fixed geometry) these 
forces vary with changes in propeller rpm, aircraft speed, induced velocity 
vector and the blade angle. How these parameters affect the forces is deter- 
mined as follows: the aircraft direction of flight and propeller rotation 

are vectorial ly combined to single velocity vector R, and the angle of attack 
ALPHA is the difference between the angular position of R and the blade an- 
gle. See Figure 2. The product of ALPHA times the slope of curves of lift 
and drag coefficients vs ALPHA is lift and drag coefficients. An examination 
of the lift and drag force equations show that these forces vary as the 
square of the velocity resultant R and directly with the lift and drag coef- 
ficients. 

These aerodynamic loads are actually pressure distributions around the face 
and camber sides of the airfoil and along the axis of the blade but they may 
be depicted as concentrated forces, the lift and drag, acting at centers of 
pressure on the airfoil cross sections and along the blade axis. The center 
of pressure on the airfoil cross section is approximately one quarter chord 
from the leading edge. The offset of this point from the pitch change axis 
of the blade retention causes aerodynamic twisting moments, ATM, about the 
blade axis. 

The center of pressure along the blade axis Is in the tip region of the blade 
and the composite lift and drag forces there produce shear forces and moments 
that are reacted at the blade retention on the hub. The effect of the shear 
forces on the hub is minute compared to the moment effects. 

If the aerodynamic loads on the blade are constant as the blade rotates, then 
the loads they produce at the blade retention on the hub are the steady bend- 
ing moment and steady ATM. 

When the aircraft climbs the aerodynamic loads on the blade become cyclic. 

The airflow into the propeller is no longer perpendicular to the plane of ro- 
tation. As the blade rotates its attitude relative to the fixed angular in- 
flow changes, that is, its leading edge at one point would be up relative to 
the flow while at 180 degrees from that point it is down. This causes 
changes in the resultant velocity vector R as well as the angle of attack as 
the blade rotates. This change is sinusoidal on each blade and the retention 
sees this as a vibratory shear and bending moment. The blade moment loads 
add at the centerline of the Prop-Fan and produce a steady moment vector and 
shear force through which the tailshaft rotates. See Figure 3. 
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FIGURE 1. HUB AND BLADE RETENTION 
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FIGURE 2. ANGLE OF ATTACK - ALPHA 
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FIGURE 3. LOADING - AERODYNAMIC-CYCLIC 




The major inertia load on the blade is the centrifugal load. It varies with 
the weight, mass center location and the square of the rpm of the blade. The 
centrifugal load acts along the blade axis and is considered a steady type 
loading. 

As the air loads bend the blade, the mass center is deflected out of the ro- 
tational plane offsetting the centrifugal load causing a restoring moment 
which must be combined with the steady aerodynamic moment. The mass center 
can be deliberately be offset to alter the effective steady moment acting on 
the retention. 

As the blade changes pitch out of the plane of rotation, part of the blade 
lies ahead of the plane and part behind. The centrifugal load acting on 
these two segments has a component of its load acting parallel but offset 
from the rotation plane. These two offset forces ahead and behind the plane 
produce a couple tending to twist the blade back into the rotation plane. 

This is the centrifugal twisting moment. This couple depends on the same 
factors as the centrifugal force. In addition the CTM varies with the sine 
of twice the blade angle. The CTM is considered to be a steady type load- 
ing. It is combined with the ATM to load the pitch change actuator system. 

The pitch change trunnion roller on the butt end of the blade is a large off- 
set distance from the retention bearing center. The twisting moment reac- 
tions on the trunnion roller produce offset moments which are included in the 
resultant steady bending moments. 

A finite element model was used to analyze the blade. The aerodynamic pres- 
sures were applied as discrete loads at the nodes of the elements comprising 
the airfoil sections of the blade. The inertia loads were applied by rotat- 
ing the mass model of the blade at the climb rotational speed about the 
Prop-Fan axis. The model was constrained at its inboard end. The blade re- 
tention and hub loads were derived from the constraint loads as forces and 
moments about an x, y, and z axis as shown in Figure 4. The loads were re- 
solved in the plane of the blade retention ball bearing centers. 

The aircraft climb condition has the highest flight vibratory loading and ac- 
cumulates a large number of cycles over the life of the aircraft. The stress 
cycles accumulate at the rate of the Prop-Fan rotational speed. The climb 
condition therefore, was used for the design case and the retention loadings 
for the climb case are tabulated in Table I. 

TABLE I. RETENTION LOADS 



Fx 

Fy 

Fz 

My 

Mz 


Newtons 

Newtons 

Newtons 

Newton-meters 

Newton-meter 


(lbs) 

(lbs) 

(lbs) 

(in-lbs) 

(in-lbs) 

STDY + VIB 

-368690 

-35240 

-14360 

-3898 

-2248 


(-82890) 

(-7922) 

(-3229) 

(-34499) 

(-19893) 

STDY - VIB 

-368690 

-29450 

-2913 

1130 

-6449 


(-82890) 

(-6620) 

(-655) 

(10003) 

(-57073) 
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The loads applied to the hub are identical to those applied to the reten- 
tion. (Figure 5 shows a hub with its associated loads). However, in the 
past, it has been found that the side loads on the retention, notably the 
torque and thrust, have a negligible effect on the hub stressing. Therefore, 
the side loads are neglected in this analysis. The loads used in the hub 
analysis come directly from Table I and are tabulated in Table II below. 
Angular values refer to number of degrees clockwise from the direction of ro- 
tation looking tip to hub. 


TABLE II. HUB LOADS 

CENTRIFUGAL FORCE 368690 newtons 

(82890 lbs) 

STEADY MOMENT (SBM) 4563 newton-meters 0 252 deg. 

(40385 in-lbs 0 252 deg.) 

VIBRATORY MOMENT (VBM) 3104 newton-meters 0 144 deg. 

(27470 in-lbs 0 144 deg.) 
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TAILSHAFT LOADS 


The loads which are applied to the tail shaft can be broken down into five 
categories: thrust, torque, engine shaft/tai 1 shaft preload, blade centri- 

fugal load, and shaft bending moment. Because of its high impressed moment 
and large number of cycles, the aircraft climb condition once again was de- 
termined to be the design limiting case. 

The thrust generated by the propeller and torque applied to the propeller are 
fed through the tail shaft. Therefore, the tail shaft must be able to with- 
stand the 32,995 newtons (7,418 lbs) of thrust and the 25,165 newton-meters 
(222,700 in-lbs) of torque the engine produces. Only a small portion of the 
blade centrifugal load is transmitted to the tailshaft, the majority being 
absorbed by the hub. The engine shaft/tai 1 shaft interface is made up of two 
cones seated in the hub and preloaded against each other by the propeller re- 
taining nut through the engine shaft. This nut is torqued 3,390 newton- 
meters (30,000 in-lbs) producing a 298,906 newton (67,200 lb) axial preload 
on the cones. The moment load applied to the tailshaft is described in 
detail in the retention load section. This is a steady bending moment that 
produces once per revolution fully reversed bending stresses in the tailshaft 
because of its rotation. The moment magnitude is 8,645 newton meters (76,500 
in-lbs). 

All of the loads, excepting the engine shaft preload, are applied where the 
tailshaft and hub are joined. The preload is applied at the cone seats where 
the engine shaft and tailshaft interface. 
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BLADE RETENTION STRESS 


The retention stresses were analyzed using a computerized analysis called 
H380. This analysis does a load balance on the bearing accounting for such 
parameters as number of balls, ball diameter, pitch diameter, forces, mo- 
ments, material properties, and geometry. The results from this analysis 
that are used in retention design include Hertzian deflection and stresses, 
position of the ball contact pattern, and moment spring rate as illustrated 
in Figure 6. For the loads listed in Table I, the following retention design 
goals have been achieved: unlimited life (>10e8 cycles) for high cycle fa- 
tigue (HCF) and a low cycle fatigue (LCF) life of 10,000 cycles as illus- 
trated in the normalized Goodman diagram in Figure 7. The predicted low cy- 
cle fatigue life of 10,000 cycles is less than the desired 50,000 cycles. 

The 10,000 cycle life is equivalent to 10,000 excursions from rest to 100% 
rpm or to 2500 pitch change cycles. This application will probably see less 
than 500 pitch change cycles or on-off cycles. This gives a five to one 
pitch change life factor or a twenty to one on-off life factor. The ball 
patterns are on the race for all flight condi ti<?ns. 

In addition to the Hertzian stress calculations for the ball race interface, 
the subsurface shear stresses were also calculated using another computerized 
analysis called P248. This calculation insures that the hardening depth ex- 
ceeds the expected peak of the subsurface shear. Figure 8 shows the depth 
characteristics of the carburizing process, and in Figure 9, the subsurface 
shear steady and cyclic stresses are plotted verses depth. The specified 
depth of .104 cm (.041 inches) for the fully hardened zone 59 Rockwell C 
minimum is sufficient to envelope the peaks in Figure 9. Figure 10 shows the 
crossover between the core material 34 Rockwell C allowables and the subsur- 
face shear which indicates a minimum depth of hardening of .297 cm (.117 
in.). The carburizing depth is defined as the point where the hardness is 
equal to the core value of 34 Rockwell C. The 50 Rockwell C equivalent depth 
is specified as .103 cm (.08 in.) to .254 cm (.1 in.). This is 2/3 of the 
full depth or .3 cm (.12 in.) min, and will J nsure that the hardness does not 
decrease to the core value above the minimum depth of .297 cm (.117 in.). 

As a further criterion of the retention design, the centrifugal load, using 
H380, was increased to 150 and 200 percent of the nominal. At 150 percent 
overspeed (the stress is plotted on Figure 7 as 150%cf), the retention balls 
and races were found to have no additional permanent deformation than the ac- 
ceptable value of .00005 in. /in. of element diameter. At 200 percent over- 
speed the force on the balls was enough to cause greater permanent deforma- 
tion. This deformation, however, will not cause a catastrophic failure in 
the retention system and meets the overspeed design requirements. 
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FIGURE 7. BLADE BEARING STRESS CRITERIA 
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FIGURE 9. SUBSURFACE RACE STRESS 
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FIGURE 10. RACE AND BALL STRESS CRITERIA SUBSURFACE OCTAHEDRAL 
SHEAR STRESS 
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BARREL STRESS 
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The stress analysis of the barrel is done by splitting it into three sec- 
tions; the front ring, the rear ring, and the barrel arms. The front and 
rear rings are analyzed separately using a computerized analysis called 
Ring-28. For this eight bladed propeller the barrel was split into six equal 
segments over the 22.5 degrees necessary to provide mirror symmetry (See 
Figure 11). Each segment is a planar cross section cut from the center line 
of the propeller radially outward. The center of gravity, area, and moment 
of inertia for each slice is computed for the front and rear rings (See 
Figure 12 and 13). These cross sectional properties, which vary with their 
azimuthal position, are then used to construct the ring. The loads, which 
originate in the barrel arms, are transposed to radial and tangential forces 
at the sections on the ring using shear flow and Wise coefficients'. The 
stresses can then be calculated for the inner and outer surfaces of the front 
and rear rings. 

The barrel arms are analyzed using a shell of revolution computer analysis 
called H727. Since the cross section must be constant, to be conservative, 
the thinnest section of the arm is used in the model (See Figure 14). The 
steady and vibratory bending moments, input as equivalent axial loads, and 
the centrifugal force, an axial force, are applied at the hub raceway to ob- 
tain the effect secondary bending through that section. 

For the top of the barrel arm, the hoop stress from Ring-28 is added to the 
hoop stress for the same position from H727 to obtain an effective combina- 
tion of the two loading patterns at that point (See Figure 11). Another re- 
latively high stressed point from the analysis is just below the ball bear- 
ing. These two worst stress points are plotted on the Goodman diagram in 
Figure 15. The amount that the point is below the line on the Goodman dia- 
gram indicates how much extra strength the part has. 

Barrel stressing was also analyzed using 3D finite element analysis, the 
results of which are summarized in Appendix A. 
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AREA - 3.51 139 
XBAR * -3.55566 
YBAR » 8.24665 
I Y-Y * 0.387344 
I X-Y * 0.463913 
I X-X * 4.30461 


AREA - 6.37453 
XBAR * 4.89617 
YBAR - 5.37813 
I y-Y * 8.80057 
I X-Y = -17.8393 
I X-X » 42.8773 


RAO SECT 00 DG 





AREA - 3.64058 
XBAR * -3.51002 
YBAR » 8.26361 
I Y-Y » 0.458043 
I X-Y * 0.618058 
1 X-X * 4.50827 




RAO SECT. 4J OG 


FIGURE 12. RING SECTION PROPERTIES 
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AREA * 4.87366 
XBAR * -3.07928 
YBAR * 8.51836 
I Y-Y = 2.02982 
1 X-Y - 2.95154 
1 X-X * 6.84723 


R Y-Y » 0.645359 
R X-X » 1.18530 
ZETA « 64.6134 
IP P-T a 8.24741 
IP X-X « 0.629647 
RP y-Y a 1.30086 
RP X-X > 0.359435 


AREA a 7.47935 > 
XBAR * 4.42450 
YBAR * 6.00063 
I Y-Y a 17.1 181 
I X-Y - -30.9509 
I X-X > 60.5261 


HUB-4 RAO SECT 13.5 PG 


XBAR « -2.37736 
YBAR - 8.90374 
I Y-Y » 9.09621 
I X-Y a 8.54895 
I X-X a 10.39558 


AREA - 9.52786 
XBAR « 3.62946 
YBAR a 6.82405 
I Y-Y - 38.1544 
I X-Y a -34.9157 
1 X-X * 85.1361 


AREA * 7.67286 
XBAR * -2.20683 
YBAR * 8.77504 
I Y-Y * 1 1.60246 
I X-Y » 7.41534 
I X-X > 7.58092 


AREA - 10.00039 
XBAR « 3.44192 
YBAR a 6.81 699 
1 y-Y » 44.4651 
I X-Y » -36.1792 
I X-X > 80.0289 


RAD SECT 18 DG 


RAO SECT 22.5 DG 


FIGURE 13. RING SECTION.PROPERTIES 
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HUB ARM 

FIGURE 14. BLADE RACE. HUB ARM MODEL 
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NORMALIZED CYCLIC STRESS 



NORMALIZED MEAN STRESS - % 

FIGURE 15. RING AND ARM STRESS CRITERIA MACHINED SURFACES 
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TAILSHAFT 


The tafl shaft was analyzed using a shell of revolution computer analysis cal- 
led H746. A radial deflection, due to centrifugal load of .035 cm. (.0014 
in.) was calculated for the vicinity of the tail shaft. The tail shaft model 
was then predeflected by this amount to reproduce the effect of its share of 
the centrifugal load. The remaining loads, as explained earlier in the tail- 
shaft load section, were applied to the nodes of the model and the stresses 
calculated are shown in Figure 16 and on the Goodman diagram shown in Figure 
17. 
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FIGURE 16. TAILSHAFT MODEL 
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NORMALIZED CYCLIC STRESS 
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FIGURE 17. TAILSHAFT STRESS CRITERIA 
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SPRING RATE 


The stress levels in the previous sections are well below allowables indicat- 
ing that there is adequate stress margin. Obviously stress was not the siz- 
ing criterion for the hub assembly. Hub wall thickness is dictated by a re- 
quired retention moment spring rate or stiffness. The blade, because of the 
more rigorous geometric constraints put on it, has a spring rate and mass 
distribution of its own which usually cannot be varied to any extent without 
upsetting the delicate balance between stresses and loads. Therefore, one of 
the major functions of the hub assembly is to provide a retention stiffness 
that will place the natural frequencies of the propeller outside the re- 
stricted range (shaded regions of Figure 18). 

The first, second, and third modes shown on Figure 18 are the first three 
blade natural frequencies. The hub in-plane and out-of-plane spring rates 
(nominally 13.7 and 15.5 million in. -lbs respectively) are plotted on these 
mode lines. The natural frequencies of the propeller system for these points 
are outside the various indicated excitation frequencies and their tolerances 
(lOt on two excitations per revolution - 2p, 7.5X on three excitations per 
revolution - 3p, and 5X on four excitations per revolution - 4p). Therefore, 
the hub in-plane and out-of-plane spring rates are acceptable for this pro- 
peller assembly. 

The spring rate of the hub assembly is divided into components which are 
evaluated separately for ease of analysis. These are: the blade shank, 
blade/race interface, race, ball, barrel arm, barrel bridge (affects inplane 
only), and barrel rings. These series springs are summed reciprocally to 
produce a total hub spring rate. 

The finite element shell of revolution model of the blade shank, race, and 
barrel arm sections provided both the stress and the stiffness values (See 
Figure 14) by computer analysis called H727. The rotation in each of the 
above mentioned parts is found for a given moment load and a spring rate is 
calculated by dividing the moment by the rotation. The ball, modeled as a 
rigid link in this analysis, to serve as part of the load path, has a non- 
linear spring rate whose stiffness is calculated by computer program called 
H380. 

The bridge consists of the hub material in between the barrel arms and thus 
its spring is effective only in the Inplane direction. The hub, as shown in 
Figure 19, is cut by planes passing through the y axis between 45 and 135 de- 
grees from the forward direction by a three dimensional computer aided design 
program called CATIA. The section properties are calculated for each slice 
(Figure 20 and 21). A graph of these section properties vs. distance along 
the bridge Is constructed (Figure 22), and a Nastran finite element model 
using beam elements is formed. A simulated moment is applied to the model, 
rotations are found and the spring rate is calculated. 

Barrel stiffness was also analyzed using 3D finite element analysis methods, 
the results of which are summarized in Appendix A. 
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FREQUENCY, HZ 



SPRING RATE, 10® IN-LB/RAD 


FIGURE 18. PROPELLER FREQUENCIES VS. TOTAL HUB SPRING RATE 








AREA * 2.81930 
XBAR * -2.90965 
YBAR * 9.16975 
£ Y-Y * 0.407504 
I X-Y * 0.158783 
I X-X > 1.46342 

HUB SECT 45 DG 



AREA » 2.33717 
XBAR > 9.27361 
YBAR - 3.08569 
I Y-Y « 1.12204 
I x-Y * -0.073189 
I X-X » 0.249207 

HUB SECT 54 DG 



AREA = 1.6771 5 
XBAR > -3.31605 
YBAR > 9.39278 
1 Y-Y = 0.098251 
I X-Y a 0.002335 
I X-X a 0.783310 

HUB SECT 72 DG 


FIGURE 20. BRIDGE SECTION PROPERTIES 
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AREA » 2.23470 
XBAR * -3.09031 
YBAR > 9.32103 
I Y-Y » 0.239029 
I X-Y a 0.088684 
I X-X - 1.001 1 1 

HUB SECT 126 DG 



AREA a 1.67715 
XBAR a .3.31603 
YBAR a 9.39256 
I Y-Y • 0.098297 
I X-Y » 0.001251 
I X-X a 0.789078 

HUB SECT 108 DG 



AREA > 2.56543 
XBAR > -2.91533 
YBAR a 9.26060 
I Y-Y a 0.370034 
I X-Y a 0.192150 
I X-X a 1.17827 

HUB SECT 135 DG 


FIGURE 21. BRIDGE SECTION PROPERTIES 



4.0 3.0 2.0 t.O 

FRONT 


0 - 1.0 - 2.0 
IN. AFT 


-3.0 -4.0 


FIGURE 22. HUB ARM BRIDGE SECTION PROPERTIES 
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The blade/race interface also has a spring associated with it. These two 
parts are loaded against each other by the pressure of the bearing causing 
Hertzian deflection between them. Approximating this condition by a cylinder 
in a groove, rotations can be calculated for the applied moments and a spring 
rate can be found. 

The barrel rings are another spring that goes into making up the total hub 
spring rate. Their deflection which make up the rotation due to applied mo- 
ments is part of the analysis done to calculate the stresses by computer 
analysis called Ring-28. The rings are loaded differently by inplane and out 
of plane moments, so, the spring rates are slightly different in the two di- 
rections . 

The last spring to be covered is the ball bearing which can readjust its con- 
tact angle under load and exhibit relatively large Hertzian deflections. 

This behavior makes its spring rate the smallest in magnitude and thus the 
most important. Because the individual spring rates are added reciprocally, 
the smallest sets the maximum limit for the total spring rate and as the 
larger Individual springs are added to it the total hub spring rate becomes 
smaller. By adjusting the bearing geometry, the ball spring rate can be var- 
ied quite substantially to obtain an acceptable total spring rate. Figure 23 
shows a survey of retentions, analyzed on our H380 bearing analysis, that 
were considered for this application. Off race indicates the ball race con- 
tact pattern edge is off the race at some loading condition and indicates un- 
suitability due to a reduction in capacity. The bearing chosen has twenty 
nine 1.588 cm. (.625 in.) diameter balls in a 15.392 cm. (6.06 in.) pitch 
diameter. 

Figure 24 shows a tabular listing of all the spring rate values and also the 
total hub assembly spring rate value. This spring rate value can then be ad- 
ded to the blade model and the natural frequencies of the whole propeller 
system can be calculated for the different blade modes. These values are 
plotted for the inplane and out of plane directions in Figure 18. This shows 
the predicted natural frequencies to be out of the range of the anticipated 
exciting forces. 
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IAMETER - INCHES 



SPRING RATE - I OES (IN-LBS)/RADIAN 

FIGURE 23. BEARING O.D. VS. SPRING RATE 
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IN PLANE 


OUT OF PLANE 


BLADE 

1*5 

1*5 

BLADE/RACE 

217 

217 

RACE 

87* 

876 

BALL 

2*. 1 29.2 

28.1 29.2 

BARREL ARM 

675 

875 

BARREL BRIDGE 

370 

— 

BARREL RINGS 

7! 

67 

TOTAL, K t 

13.27 14.03 

15.04 16.03 


K t = 1/(1/K,+ 1/K n ) 


FIGURE 24. 

HUB SPRING RATE IN* 

LB/RAD* 10 6 
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WEIGHTS 


The weight breakdown for the hub assembly parts appears in Figure 25. 
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POUNDS/PROP 


HUB 218.6 

RET. 

BLADE/RACE 4.48 

BLADE SEAL 2.08 

BLADE RING 16.20 

BALL (232) 8.4 

BALL SEPARATOR 0.80 

FRICTION STRIP 0.16 


32.1 2 

TOTAL 250.8 

FIGURE 25. HUB AND RETENTION WEIGHT 
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CONCLUSION 


The hub assembly meets all design requirements. All components will meet or 
exceed life expectations. The single exception is the retention bearing low 
cycle fatigue life which is predicted to be less than 50,000 cycles but will 
exceed the anticipated utilization by a twenty to one factor. 
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APPENDIX A 

HUB RETENTION STIFFNESS AND STRESSING 
3D FINITE ELEMENT ANALYSIS 
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Introduction 


A 3D finite element analysis of the LAP barrel was performed to verify the 
retention stiffness and barrel stressing for centrifugal, 2P and IP loading 
conditions. Solution accuracy checks on model mesh size and constraints were 
also made. The following sections summarize the FEA results and compares the 
FEA results with the ring analysis. 

Concl usions 

1. The FEA in-plane (IP) barrel stiffness is 37 x 10 6 in-lb/rad. The 
out-of-plane (OOP) barrel stiffness is 70 x 10 6 in-lb/rad. 

2. The FEA IP stiffness is 33% lower than the stiffness obtained from the 
ring analysis. The FEA OOP stiffness is 15% higher than the ring 
analysis. 

3. Stress levels obtained for all the FEA load cases are below the 
allowables. 

4. The barrel stiffness is not sensitive to the model ball load application 
point. 

5. For the 2P moments of 10,000, 20,000 and 28,000 in-lbs. the IP retention 
stiffness decreased a maximum of 3%. The OOP stiffness increased only 
.3%. 

6. In this analysis, the 8-noded brick finite element model is adequate for 
determining barrel stiffness and stressing. 
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Methodology 


F. E. Model : Structural symmetry permitted the modeling of only 

one-sixteenth of the barrel; see Figure A1 . The model consists of 251 
8-noded <8N) CHEXA and six-noded (6N) PENTA solid elements. 

The mesh density around the barrel arm was defined so that a race node would 
represent a ball center. Nominal dimensions were used to construct the model. 

To prevent rigid body motion, the end of the barrel tailshaft, which butts up 
against the thrust face of the engine, was rigidly fixed. The boundary 
conditions on the symmetry faces of the barrel are automatically set in the 
cyclic symmetry solution sequence to reflect the applied load condition. 

Loading 


Sinusoidal ball load and bearing race contact angle distributions were 
obtained for each load case using the Nonlinear Single Race Retention 
Stiffness Analysis Program F264CL and a Fortran program. Unlike current 3D 
FEA barrel analyses in which the ball loads were distributed over several 
nodes at each ball center location, the ball loads in this analysis were 
applied to a single race node at each ball center location; see Figure A2. 
This method of loading recognizes that the mesh density in the race area is 
too coarse to approximate the Hertzian ball load distribution and that the 
local "deflection" footprinting effect is negligible. This was verified by 
running a case with two loaded nodes. Asymmetric arm loadings for 2P and IP 
modes are accounted for by using the Dihedral Cyclic Symmetry solution format 
in MSC/NASTRAN. 

Stiffness 


In- and out-of-plane retention stiffness values are obtained by performing a 
series of applied load iterations on the barrel. These iterations are 
necessary due to the nonlinear action of the barrel under load. The 
iteration process starts by initially applying a uniform load distribution 
with a constant contact angle to the raceway. Axial and radial displacements 
are obtained at the top, bottom, in-plane left and right locations on the 
race at the appropriate contact angle; see Figure A3. Compliances are then 
calculated at these points. These compliances are added to the blade, 
blade/race, race and ball without the kinematics, are input into F264CL to 
represent the entire retention flexibility. The kinematics are included in 
F264CL. This results in a new load distribution which is then applied to the 
FEA model. The process continues until two consecutive Iterations result in 
compliances which are within 5% of each other. In- and out-of-plane cocking 
angles are then obtained from F264CL and converted into retention moment 
stiffnesses. 
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FIGURE At. 3-D FINITE ELEMENT MODEL - 1/16 SEGMENT 






POINTS A, B, C AND D ARE 
DISPLACEMENT RECOVERY POINTS 





FIGURE A3. DATA RECOVERY POINTS FOR STIFFNESS ITERATION PROCESS 
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Analysi s 


Stiffness Results : A centrifugal force of 86,000 lbs. was applied to the 

model. This load application resulted in an in-plane retention stiffness 
<IP) of 10 x 10 6 in-lb/rad. and an out-of-plane retention (OOP) stiffness 
of 18 x 10 6 in-lb/rad. From these values in- and out-of-plane barrel 
stiffness terms are obtained and further reduced down into their fundamental 
components of the ring, arm and bridge to determine relative component 
stiffnesses; see Figure A4. The bridge is the softest spring in the barrel. 
The magnitude of these stiffnesses, however, are approximate in that 
assumptions are made where the arm ends and the ring begins; see Figure A5. 

A comparison between the FEA and the ring analysis is shown in Figure A6. 
Differences between the two approaches are primarily due to the use of a 
sinusoidal load and contact angle distribution in the finite element analysis. 

The 2P bending retention stiffnesses were obtained for a centrifugal force of 
86,000 lbs. plus 10,000, 20,000 and 28,000 in-lb. vibratory moments. For 
each load case, the applied moment was split into IP and OOP components where 
the resultant moment was oriented at the calculated IP bending mode 
orientation of 0 - 142°; see Figure A7. This orientation represents a 
flatwise bending on the blade. A summary of the IP and OOP bending retention 
stiffnesses are shown in Figure A8. An arm load phasing of .707 on each arm 
was found to produce the lowest barrel stiffness and highest stresses. For 
the range of applied moments, the IP and OOP stiffnesses changed a maximum of 
3% and .3% respectively. Based on this result, it is assumed that inclusion 
of the steady bending moment would not significantly reduce the retention 
stiffnesses . 

In- and out-of-plane retention stiffnesses for the centrifugal force and 2P 
cases are plotted in Figure A9 for the purpose of obtaining flatwise and 
edgewise retention stiffness values for the 2P case. The edgewise and 
flatwise retention stiffnesses on the centrifugal force only line were 
obtained by calculating barrel compliances at points A, B, C and D which lie 
on the edgewi se-flatwi se axes and following the load and stiffness procedure 
previously described. A surve was fit to the four points. A similar curve 
was assumed for the 2P data. A flatwise retention stiffnesses of 1 1 x 10 6 
in-lb/rad. and an edgewise retention stiffness of 13.5 x 10 6 in-lb/rad. are 
obtained from this procedure. 

A IP flatwise loading retention stiffness was also obtained for a centrifugal 
force of 86,000 in-lb. and a 28,000 in-lb. vibratory moment at 0 * 142°. 

The load arm phasing is shown in Figure A10. This loading resulted in an OOP 
stiffness of 5 x 10* in-lb/rad. and an IP stiffness of 6 x 10 6 in-lb/rad. 
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COMPONENT 

OOP (IN-LB/RAD) 

IP (IN-LB/RAD) 

RING 

172 X 10 6 

191 X 1 0 6 

BRIDGE 

- 

46 X 10 6 

ARM 

118 X 10 6 

2 X 10 n 

TOTAL 

70 X 1 0 6 

37X10 6 


FIGURE A4. BARREL MOMENT STIFFNESSES 
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COMPONENT 

OOP (IN-LB/RAD) 

IP (IN-LB/RAD) 


A 

B 

A 

B 

RING 

172X10* 

95 X 10® 

191 X 10 6 

97 X 10® 

BRIDGE 

— 

— 

46X10* 

60 X 10® 

ARM 

118X10* 

267 X 10* 

2X10*1 

4 X 10*0 

TOTAL 

70 X 10 6 

37 X 10 s 


FIGURE A5. EFFECT OF POINT SELECTION FOR BARREL COMPONENT 
STIFFNESS DETERMINATION 
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FEA 


COMPONENT 

OOP (IN-LB/RAD) 

IP (IN-LB/RAD) 

RING 

172 X 10® 

191 X 10 6 

BRIDGE 

— 

46 X !0 6 

ARM 

11 8 X 1 0 6 

2 X 10 1 1 

TOTAL 

70 X 10® 

37X10 6 


DESIGN 


COMPONENT 

OOP (IN-LB/RAD) 

IP (IN-LB/RAD) 

RING 

67 X 1 0® 

71 X 1 0 8 

BRIDGE 

— 

370 X 10 6 

ARM 

675 X 10« 

675 X 10 6 

TOTAL 

61 X 10 6 

55 X 1 0 6 


FIGURE A6. COMPARISON OF FEA AND DESIGN BARREL STIFFNESS RESULTS 















FIGURE A7. 2P APPLIED LOADS 


58 



• 2P - PHASING - 0.707 



APPLIED LOADS 

OOP (IN-LB/RAD) 

IP (IN-LB/RAD) 

C.F. ± M 

86,000 + 10,000 

16.87 X 10® 

7.94 X 10® 

86,000 + 20,000 

16.93 X 10® 

7.80 X 10® 

86,000 + 28,000 

16.92 X 10® 

7.69 X 10® 


FIGURE A8. IN-AND OUT-OF-PLANE BENDING 

RETENTION STIFFNESS SUMMARY FOR 2P MODE 
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0° 38° 32° 

90° 

180° 

360° 

OOP (TOP) 

IF 

OOP (BOTTOM) 




(0) AZIMUTHAL ARM LOCATION 



FIGURE A9. 2P FLATWISE AND EDGEWISE RETENTION STIFFNESS DETERMINATION 
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ARM 1 


8 



DEFLECTED SHAPE 
IS SUPERIMPOSED 



0 


142 ° 



ARM S 


FIGURE A10. IP LOAD PHASING 


61 


Stress Results 


Average nodal point equivalent (Von Mises) stresses were evaluated for the 
centrifugal force, 2P and IP load cases. Stresses in these cases are 
predominantly due to the centrifugal force. Stress contours for this case 
are shown in Figure All. Peak stress contours are highlighted. 

Cyclic stresses resulting from the IP load case are shown in Figures A12 and 
A13. Both the right and left segments of arms 1 and 5 are shown because of 
the asymmetric applied loading. Due to the nonsymmetric response of the 
barrel, the cyclic stress does not fully reverse. In all the cases surveyed, 
the stress levels are below the allowables for steel. 

Analysis Checks 

Several checks were made on the analysis procedure to determine solution 
sensitivity. The centrifugal force load case was used for these checks 
because it required the least amount of solution and evaluation time. The 
first check made was with, respect to the selection of the load point on the 
race. This was done to determine the sensitivity between the load point and 
the resultant stiffness. Figure A14 shows the two load point locations 
selected. The first choice is the location used throughout the study. It 
was felt that this race location was a fair representation of where the ball 
bearing would be located during the various loadings. The second choice was 
assumed to be an extreme case that could not occur under normal operating 
conditions. The load points selected resulted in stiffness differences of 
less than 21 indicating a low sensitivity to the load point selection. 

A solution convergence check was made by constructing a model using 20-noded 
bricks Instead of 8-noded bricks. This Increased the number of degrees of 
freedom from 1,300 to 4,600 without changing the number or distribution of 
the bricks. A comparison of displacements between the 8-noded and 20-noded 
model for a centrifugal loading Is shown in Figure A15. The 20-noded model 
results in higher displacements. This translates into an out-of-plane 
stiffness which is approximately 2% lower than the 8-noded model and 61 lower 
in-plane. Comparison of stresses shows, as expected, a greater difference; 
see Figure A16. However, due to the placement of the midside nodes in the 
20-noded brick, errors exist in the stresses. Since the stress levels are 
generally low with regard to the allowables, the stress results do not 
present a problem. 

Summary 

The use of 3D finite element analysis coupled with the program F264CL has 
provided a significant amount of insight into the structural action of the 
barrel. Both the stiffness and stress levels are within the allowables. 
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FIGURE All. CENTRIFUGAL FORCE EQUIVALENT STRESS 





RIGHT SEGMENT (MAX LOAD) 


28,000 IN-LB 4 52° OFF OOP 



65 


FIGURE A1 3. IP - CYCLIC STRESSES ON FULL MOMENT ARM 
LEFT SEGMENT (MIN LOAD) 


11.46 X 10 6 16.84 X 10 6 


IWWr 


11.70 X 10* ‘1 6.92 X 10* 


FIGURE A1 4. STIFFNESS SENSITIVITY TO APPLIED LOAD POINT 
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FIGURE A15. COMPARtSON OF 8-NODE D VERSUS 20-NODED MODELS 




FIGURE A16. COMPARISON OF 8-NODED VERSUS 20-NODED MODELS 



APPENDIX B 

HUB ASSEMBLY LAYOUT DRAWING L-1 4325-5 
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